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ABSTRACT 

We report the clear detection of the 2175 A dust absorption feature in the optical afterglow spectrum 
of the gamma-ray burst (GRB) GRB 070802 at a redshift of z = 2.45. This is the highest redshift 
for a detected 2175 A dust bump to date, and it is the first clear detection of the 2175 A bump 
in a GRB host galaxy, while several tens of optical afterglow spectra without the bump have been 
recorded in the past decade. The derived extinction curve gives Ay = 0.8-1.5 depending on the 
assumed intrinsic slope. Of the three local extinction laws, an LMC type extinction gives the best 
fit to the extinction curve of the host of GRB 070802. Besides the 2175 A bump we find that the 
spectrum of GRB 070802 is characterized by unusually strong low-ionization metal lines and possibly 
a high metallicity for a GRB sightline ([Si/H] = -0.46 ± 0.38, [Zn/H] = -0.50 ± 0.68). In particular, 
the spectrum of GRB 070802 is unique for a GRB spectrum in that it shows clear C I absorption 
features, leading us to propose a correlation between the presence of the bump and C I. The gas to 
dust ratio for the host galaxy is found to be significantly lower than that of other GRB hosts with 
N(H I)/ Ay = (2.4 ± 1.0) X 10^^ cm"2 mag"\ which lies betwee n typical MW and LM C values. Our 
results are in agreement with the tentative conclusion reached by 'Gordon et ajj ()2003[ ) that the shape 
of the extinction curve, in particular the presence of the bump, is affected by the UV flux density in 
the environment of the dust. 

Subject headings: dust, extinction - galaxies: ISM - gamma rays: bursts - galaxies: abundances — 
galaxies: distances and redshifts 



1. INTRODUCTION 

Dust extinction curves quantify as a function of wave- 
length the amount of light 'lost' due to scattering and 
absorption of the light by dust particles along the line 
of sight from an object to the observer. In the Milky 
Way (MW), the extinction curve has been extensively 
mapped and has been shown to follow an empirical 
single-parameter fun ction for almost all lines of sight 
(jCardelli et al.|[T"98i/ ). The most characteristic feature of 
this function is a broad 'bump' (i.e., excess extinction) 
centered at 2175 A first discovered bv IStecherl (|1965f ). 
Today, more than 40 years after its discovery, the origin 
of the feature remains unknown although several candi- 
dates have been suggested (see i) ll.2p . 

As outlined below (in § ll.3|) , a promising method for 
studying extinction curves in distant galaxies is to use the 
spectra of the afterglows of gamma-ray bursts (GRBs) 
as the backlight against which the extinction curve can 
be inferred. The intrinsic spectral energy distributions 
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(SEDs) of GRB afterglows are very simple, consisting of 
a single or a broken power-law from the X-ray band to 
the infrared (IR) and hence it is relatively easy to in- 
fer the shape of the extinction curve by measuring the 
curvature of the optical spectrum, the deviation of the 
afterglow spectrum extending from the X-ray extrapola- 
tion, or from the X-ray-IR interpolation. 

In this paper we present the detection of the 
2175 A bump in a GRB absorption sys tem at z = 
2.45. The detection was briefly reported Ijv lFvnbo et al.l 
(|2007D , but this paper gives the detailed analysis of our 
detection. The detection of the bump has also sub- 
sequently been conf irmed based on photometry alone 
(|Kriihler et al.|[2"008f ). This is the first time that a clear 
signature of the 2175 A feature has been observed in a 
GRB host galaxy and the highest redshift detection of 
the bump to date. 

The paper is organized as follows: In the remainder of 
this section we discuss the nature and previous detections 
of the 2175 A feature (§ [OD and the use of GRBs to 
infer extinction curves (§ II. 3p in more detail. In § [5] we 
present the observations and data analysis. We discuss 
the properties of the derived extinction curve in § [3l The 
results and possible tracers of the bump are presented in 
§ HI and finally we summarize our conclusions in § [Sj 

1.1. Extragalactic extinction curves 

The study of extinction curves outside the Milky Way 
is a challenging task and has mostly been limited to 
galaxies in the Local Group. While the M31 extinction 
curve resembles that of the Milky Way (|Bianchi et al.l 
Il996l ). the extinction curves of the Large and Small 
Magellanic Clouds (LMC and SMC) show significant 
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variation. Lines of sight in the LMC can be broadly 
put into two categories, one having extinction curves 
similar to that of the Milky Way ("LMC average") 
and the other showing extinction curves with a much 
less prominent bump and a steeper rise into the UV 
/"LM C2 supershell" iNandv et all [ToSlI : iMisselt et all 
fl999tlGordon et al.ll2003D . While the phrase "LMC type 
of extinction" refers to extinction curves of the second 
type (i.e. less prominent bump and a steep UV rise), it 
should be kept in mind that Milky Way sightlines also 
exist in the LMC. Similarly, for the SMC, the canoni- 
cal SMC extinction curve shows no e vidence for a bump 
and an even s t eeper rise into the UV (jPrevot et al.lll98-i : 
iGordon et al.l I2003D although one out of five lines of 
sight shows the 2175 A bump ( Lequeux "eTallflQSl ). It 
is also worth noting that there are a few known lines 
of sight in the Milky Way with LMC type of extinction 
(jClavton et al.l [2OO0I) and a single line of sight with a 
measured SMC type of extinction (jValencic et al.l[2003h . 

Relatively little is known about dust extinction out- 
side the Local Group as the method of measuring ex- 
tinction curves by comparing the spectra of single stars 
is not applicable at larger distances. Several methods 
have been attempted, including studying statistica,! sam- 
ples of red dened and standard quasars (jPei et al.lll999l : 
IVladilo eFal. 2008). studying individual supernova (SN) 
la lightcurves and using gravitationally lensed quasars. 
A recent o verview of these s t udies and their results can 
be found in lElfasdottir et al.l (|2006D . The SEDs of GRBs 
may also be used, as discussed in § 11.31 These studies 
have found varying types of extinction which can deviate 
significantly from MW-type extinction, leaving open the 
question of how common MW type dust is. In particu- 
lar, prior to the present study there have only been three 
robust detections of the 2175 A bump in individual ex- 
tragalactic systems beyond the local group. The first is 
in a l ensing galaxy at a redshift of z = 0.83 (jMotta et al.l 
I2OO2D . Th e second is from a dain ped Lya system at 
z = 0.524 (jJunkkarinen et al.|[200^ . The third is in an 
intervening ab sorber at z = 1.11 toward GRB 060418 
(jEllison et a l.| [2006 ). 

Extinction curves are the prime diagnostic tool avail- 
able to study dust in the optical/UV regime. They 
depend sensitively on the composition of the dust 
(jHenning et al.l [20041 ). allowing estimates to be made of 
the nature and the origin of dust and its dependence on 
cosmic time, metallicity, etc. Furthermore, in studies of 
distant objects, e.g., using SNe la as standard candles or 
determining flux ratios of strongly lensed QSOs, it has 
become increasingly important to accurately determine 
the amount and wavelength dependence of the dust ex- 
tinction. 

1.2. The nature and origin of the 2175 A extinction 
feature 

iFitzpatrick fc Massal ()1986f) studied the 2175 A inter- 
stellar extinction bump in the direction of 45 reddened 
stars in the Milky Way, and found that it displays a peak 
whose central wavelength Aq is remarkably constant with 
extreme deviation of only ~ ±17 A from the mean po- 
sition of Aq — 2174.4 A. As this is significantly larger 
than the measurement uncertainty it seems to indicate 
a real, but small, variation of the peak position. The 
full width at half maximum (FWHM) of the bump has a 



larger range of intrinsic values, from 360 to 600 A. There 
seems to be no correlation between the width and the 
central wavelength of the bump. 

The fact that the 2175 A feature becomes progressively 
fainter in the MW, LMC and SMC has been attributed 
to the pro gressively lower m etal abundances of the LMC 
and SMC ^Fitzpatrickll2004D or to differences in the ra- 
diative environment "("Gordo n et"aLlll997t [Clavton et all 
l2000t iMattsson et a l. 2008V IGordon et al.l (|2003[ ) stud- 
ied the differences in the extinction of the MW, SMC and 
LMC finding tentative evidence that the bump strength 
c orrelates wi t h dus t-to-gas ratio. 

iNoU et all (|2007l ) found in their study of 108 massive, 
UV-luminous galaxies at 1 < z < 2.5 that there is a cor- 
relation between heavy reddening and the presence of a 
2175 A feature. The least reddened objects have SEDs 
consistent with the average featureless and steep extinc- 
tion curve of the SMC. For their objects at 1 < z < 1.5 
a significant UV bump is present in galaxies that ap- 
pear disk-like and which host a rather large fraction of 
intermediate-age stars (i.e., from 0.2 to 1-2 Gyr old). 
More clumpy and irregular objects seem to have extinc- 
tion curves resembling those of nearby starburst galaxies. 

Several candidates have been proposed as the dust par- 
ticles responsibl e for the bump ranging from iron-poor 
silicate grains (jSteel &: DulevI |T987) to carbon aceous 
materials such a.s carbon-onions (Henra rd et al.iri 997). 
graphite grains (jStecher fc Donni ,1965. : .Draindligs g) or 
the polycyclic aromatic hydrocarbons (PAHs) believed 
also to be responsible for the prominent broad emission 
features fou nd in the mid-infr ared ( Dulev & Scahra 1998); 
lDulevll2006l : lCecchi-Pestellini et al.li2008D . 

Graphite has been considered a very promising, though 
contentio us, candidate for explai ni ng the 2175 A bump 
(e.g. IStecher fc DonnI 11 965': 'Fitzpatrick fc Massal 
1986t iMathisI 119941: iVosh chinnikov 19 90: Sorrell 
19901: iDraine fc Malhotral [1991 



, . iRouleau et al.l 

IWill fc AannestadI 119991 : lAndersen et al.l l200a 



Clavton et al.l 



20031) . In an extensive investigation, 

Draine fc Malhotral (|1993l ) conclude that if graphite 
particles are the carriers of the 2175 A peak, a variation 
in their optical properties must be present, e.g., as a 
result of varying amounts of impurities, variations in 
crystallinity, or changes in its electronic structure due to 
surface effects. The observed lack of correlation between 
the central wavelength and the FWHM of the peak is 
therefore a challenge for the hypothesis that graphite 
particles are the source of this peak. 

Large PAH molecules are expected to at least con- 
tribute to the 2175 A feature as the interior carbon 
atoms have electronic orbitals closely resembling those 
in graphene and an oscillator strength f or each C ex - 
pected to be close t o the value for graphite (|Drainell2003D . 
iJoblin et al.l (|1992f l have suggested that PAHs could be 
entirely responsible for the feature, which has been ques- 
tioned by .Mathis (1994) due to the lack of an absorption 
feature at about 3000 A in the spectra of PAHs. Another 
argument against the PAHs as the main contributers is 
that the constancy in wavelength along with a rather 
wide variation in the width, seems surprising if the bump 
is caused by a mixture of widely varying materials. Such 
properties are more likely to occur as a result of coat- 
ing (or surface hydrogenation) of carbon grains (jMathisI 
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[ToMlDrainelFLil [20071 : ICecchi-Pestellini et al.l[200i) . 

Finally, iron-poor silicates in the form of partially hy- 
drogenated amorphous Mg2Si04 particles ha ve been sug- 
geste d as the carriers of the 2175 A peak (St eel fc Dulevl 
[1987). The absorption is observed in silicates of greatly 
differing composition and under a wide range of different 
sample prepar ations and outgasing conditions. However, 
iMathisI ([19961 ) indicated that the silicates may be prob- 
lematic as carriers. 

It seems most plausible that the carrier of 2175 A 
peak is a composition of graphitic and amorphous car- 
bon phases, PAHs and glassy iron poor silicates. Vari- 
ous models ha ve been proposed with diff erent composi- 
tions (see e.g.. lMathiill996t lDrainell2003l and references 
therein). 

1.3. GRB afterglows and extinction 

GRB afterglows have very simple intrinsic synchrotron 
spectra which are essentially featureless power laws 
{fi, oc v^^) across the wavelength range of interest for 
measuring extinction curves (e . g.. [ Galama et al.|[l998[ : 
iGranot fc Sari 200^ ISari et"an il998D. This is a signif- 



icant advantage over using QSOs or galaxies as back- 
ground objects. In cases where the GRB physics places 
the cooling frequency between the optical and X-ray 
regime (at the time of observation) there will be a cool- 
ing break in the spectrum leading to a change in spec- 
tral slope of A/3 = 0.5 (softer spectrum in the X-rays). 
The X-ray slope is usually around 1.0-1.2 (see e.g., 
[Liang etaP l2007l ) and so the intrinsic optical slope is 
typically between 0.5 and 1.2. Moreover, the optical-NIR 
range is typically well-described by a single power-law 
s egment. 

iReichar^ ([1998f ) first attempted to c onstrain the ex - 
tinction towards GRB 970508 (see also [Reicharll2001[ ). 
Several subsequent studies confirmed that a simple 
SMC-type/linear extinction model provid ed good fits 
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ings of a low dust-to-gas ratios in GRB afterglows ex- 
hibiting damped Lya absorption in the GRB host galaxy 
supports the apparent prevalence of SMC-type dust 
in GRB environments (pensen et al.|[200il : [Hiorth et al.l 
[200l iVreeswiik et al.ll2004D . 

Gray (wavelengt h -independent) ext i nction laws 
jGalama et al.' 2003; 'Strat ta et all [2007t [Perlev et al" 
2008: Nardini et al. 20(38.) or MW dust ([Kann et al 
20061 : [Schadv et a'ir[2007[ ) have been suggested 



m rare 

instances for GRB environments although it should be 
noted that these claims are somewhat model dependent 
(gray dust) or statistically insignificant (MW-type 
dust). The 2175 A feature itse lf has been searche d for 
and ruled out i n GRB 050401 ([Watson et al.[ [20061). In 
GRB 000926 ([Fvnbo et al.l [2001aD and GRB 020124 
([Hjorth et aH [20031 ) SMOtvpe extinction was preferred 
although the existence of a 2175 A bur np was not sam- 
pled b y the broad-band observations. [Vreeswiik et al.l 
i[2006f) found a bump in the spectrum of GRB 991216 
but at 2360 A if at the rcdshift of the GRB {z = 1.02). 
To be consistent with restframe 2175 A it would have 
to be due to an even more distant absorber at z = 1.19 



which was however not detected as an absorption system 
in the afterglow spectrum. The first (and thus far, only) 
clear detection of the 2175 A bump in the foreground 
of a GRB, consistent with a detect ed intervening Mg I I 
absorber (z — 1.1), was reported bv [Ellison" "eiTan ([20061 ). 

A signific ant fraction of GRB afterglows are opti- 
cally 'dark' ([Fvnbo et al.[ [2001bl: [Jakobsson et al.|[2004a[ : 
[Rol et ahl [2005f ). It is likely that a large fraction of 
these afterglows go undetected in the opt ical because 
of extinction ( e.g., Diorgqvski et al.[ [20C)1[). Recently 
[Jaunsen et all ([20081 ) and [Tanvir et all ([20081) demon- 
strated the validity of this dust obscuration hypothesis 
by detecting two highly reddened systems. These stud- 
ies highlight the difficulty in probing systems with sig- 
nificant extinction at high redshift: GRBs with high ex- 
tinction are generally not accurately localized (because 
of the lack of an optical afterglow position) and there 
are no detailed constraints on their extinction laws, par- 
ticularly in the UV, for systems with of order Ay ^ 5. 
Such selection effects must be kept in mind before mak- 
ing statistical inferences from extinction laws of GRB 
afterglows. 

Studies using the X-ray extrapolation to set limits on 
the total extinction ([Vreeswiik et al.[|l999f ) or using the 
optical spectral curvature to determine the reddening 
([Ramaprakash et al. 19981 have now been carried out for 
nearly a decade. These studies may be further improved 
if one obtains the largely extinction- free IR flux, typically 
needing observations in the mid-IR, a nd is the subject o f 
our ongoing Spitzer campaign (see also [Heng et al.l20"08[ ). 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. VLT imaging, spectroscopy and redshift of GRB 

070802 

GR B 070802 (iBarthelmv et al.l 120071) was detected by 
Swift (iGehrels et al.[[2004r i on 2007 Aug 2.29682 UT. The 
ESO VLT started observations about one hour after the 
burst in Rapid Response Mode, starting with a photo- 
metric sequence spanning from the _B-band through V, 
R, and / to the z-band in excellent seeing conditions (< 
0.5 arcsec) and at low airmass (< 1-4). We detec t ed the 
source previously detected bv [Berger fc Murphvl ([2007f ) 
in all bands (see Table [1]). The source was very red 
an d classifies as a dark bur st according to the definition 
of [Jakobsson et all ([2004a[ ). Based on 5 R-band points 
taken from about 1 to 3 hours after the trigger we in- 
fer a decay slope of 0.62± .05 co nsistent with the mea- 
surement in [Kriihler et al.[ ([2008f ) at the same time. We 
infer a celestial position of R.A.(J2000) = 02'^26"35'!77, 
decl.(J2000) = -55°31'39'.'2 (calibrated against USNO- 
Bl.O). Following the afterglow detection, spectroscopy 
was immediately started, using the F0RS2 instrument 
equipped with the grism 300V (- 10 A FWHM spec- 
tral resolution). A total of 5400 s of exposure were ac- 
quired, split into three exposures, with mean epoch of 
Aug 2.38935 UT (2.2 hours after the burst). The spectra 
were obtained in excellent seeing conditions (about 0'.'5) 
and a relatively small airmass of 1.2. We used a I'.'O slit 
and therefore slitloss should be neglicible. The spectra 
were flux calibrated using a spectrum of the spectropho- 
tomctric standard star LTT1020 obtained on the same 
night. 

The optical spectrum of the afterglow is shown in 
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TABLE 1 

Log of photometric observations 



Mean epoch 


Time since 


Exposure 


Filter 


Mag. 


(2007 Aug. UT) 


trigger (hr) 


times (s) 




(Vega) 


2.33589 


0.94 


60 


2 


20.39±0.05 


2.33753 


0.98 


40 


I 


21.36±0.10 


2.33443 


0.90 


30 


R 


21.86±0.06 


2.34353 


1.12 


30 


R 


22.03±0.06 


2.34451 


1.14 


60 


R 


22.18±0.06 


2.35186 


1.32 


60 


R 


22.26±0.05 


2.42306 


3.03 


60 


R 


22.72±0.05 


2.33881 


1.01 


40 


V 


22.69±0.15 


2.34066 


1.05 


150 


B 


23.99±0.28 



GRB070802 




4000 5000 6000 7000 8000 
Observed wavelength (A) 



Fig. 1.— The flux-calibrated spectrum of the GRB 070802 after- 
glow vs. observed wavelength. Metal lines at the host redshift are 
marked with solid lines whereas the lines from the two intervening 
systems are marked with dotted lines. The broad depression cen- 
tered around 7500 A is caused by the 2175 A extinction bump in 
the host system at 2abs = 2.4549. The cleaned spectrum used for 
the extinction curve analysis in § [3] is overplotted in red. Telluric 
features are indicated with ©. 

Fig. [TJ The afterglow spectrum is characterized by a red 
continuum with a broad dip centered at Aobs ~ 7500 A. 
There is a break in the spectrum at Aobs ~ 4250 A blue- 
ward of which little residual flux, if any, is detected. We 
identify the latter feature as the imprint of the Lya forest 
of absorption lines at Zabs < 2.5. Numerous narrow ab- 
sorption lines are detected throughout the spectrum (see 
Table [2]). The highest redshift absorber is a strong metal 
line system at Zabs = 2.4549, which we adopt as the GRB 
host galaxy redshift. In addition we detect two interven- 
ing Mg II absorption systems at z = 2.078 and z = 2.292. 
The system at z — 2.078 is very strong with a restframe 
equivalent width of the 2796 A line of Wr = 3.9 A. 

2.1.1. HI content and metal lines 

We measured the total neutral hydrogen column den- 
sity (in units of cm~^ throughout) of the absorber at 
the GRB redshift by fitting a damped Lya absorption 
line profile fixed at the redshift of the detected metal 
lines. Both the saturated core and the red damped wing 
of the Lya line are constrained by the residual flux at 
Aobs > 4250 A. We derived log7V(H I) = 21.5 ± 0.2 (see 
Fig. El). 

As listed in Table^l a large number of metal species are 
identified at the GRB redshift. The detected lines from 
singly ionized elements are very strong compared to other 
GRB- or QSO-DLA sightlines. For instance, Wr = 2.9 A 
for Si IIA1526 is larger than for any LBG/GRB/QSO- 



TABLE 2 

Absorption lines in the afterglow spectrum of 
GRB 070802. 



^obs [A] 


-'^rest iM 


z 


Feature 


Wo [A] 


4200.0 


1215.7 


2.455 


Lya 




4607.9 


1335.3 


2.4529 


C II / C II 


* 16±6 


5273.9 


1526.7 


2.4544 


Sin 


10±2 


5350.9 


1548.2/1550.8 


2.4533 


C IV 


7.3±1.8 


5390.7 


1560.3 


2.4549 


C i'' 


2.4±1.2 


5556.4 


1608.4 


2.4546 


Fe II 


10.5±1.6 


5724.5 


1656.9 


2.4550 


C I 


4.7±1.1 


5771.2 


1670.7 


2.4544 


Al II 


14.3±1.2 


6246.5 


1808.8 


2.4549 


Si II 


3.9±1.0 


6404.9 


1854.7 


2.4533 


Al III 


2.4±1.0 


7000.1 


2026.1 


2.4550 


Zn II 


4.2±1.0 


7000.1 


2026.5 


2.4553 


Mgi 


7104.1 


2056.3 


2.4548 


Cr II 


2.0±0.9 


7124.8 


2062.2 


2.4545 


Cr II 


2.53±0.9 


7126.3 


2062.7 


2.4545 


Zn II 


7810.7 


2260.8 


2.4549 


Fe II 


6.6±1.0 


8094.1 


2344.2 


2.4515 


Fe II 


16.5±2.0 


8200.5 


2374.5 


2.4536 


Fe II 


12.7±1.5 


8228.0 


o o o o o 

2382.8 


2.4531 


Fe II 


18.8±1.5 


8277.1 


2396.4 


2.4540 


Fe II* 


2.9±1.4 


8308.7 


2405.6 


2.4539 


Fe II** 


2.8±1.4 


8902.8 


2576.9 


2.4549 


Mnll 


5. Oil. 2 


8934.4 


2586.7 


2.4538 


Fe II 


19.5±1.3 


8979.5 


2600.2 


2.4534 


Fen 


29.0±1.3 


8515.5 


2586.7 


o onon 

z.zyzu 


Fe II 


2.3it0.6 


8560.6 


2600.2 


2.2923 


Fen 


2.2±0.6 


9208.6 


2796.3 


2.2931 


Mgii 


1.8±0.6 


9226.5 


2803.5 


2.2911 


Mgii 


1.4±0.6 


5141.4 


1670.7 


2.0774 


Alii 


6. Oil. 7 


7218.7 


2344.2 


2.0794 


Fen 


3.4±1.0 


7311.4 


2374.5 


2.0791 


Fen 


4.5±1.0 


7335.7 


2382.8 


2.0786 


Fen 


4.9±1.0 


8005.9 


2600.2 


2.0790 


Fen 


4.7±1.1 


8607.0 


2796.3 


2.0780 


Mgii 


11.3±0.6 


8631.0 


2803.5 


2.0787 


Mgii 


11.5±0.6 


8780.1 


2853.0 


2.0775 


Mgi 


2.7±1.1 


Note. — 


The equivalent widths are given in 


the observer 



frame, Wo- 
The line is blended with C I* and C I** 




Relative velocity (km s ^) 



Fig. 2. — A section of the spectrum centered on the position of 
Lya at z^bs = 2.4549. Overlaid is the best fitting DLA profile and 
associated 1 cr uncertainty, corresponding to log iV(H I)= 21.5±0.2. 

DLA observed to dat e, hinting at a high metallicity 
(jProchaska et al.ll2008l ). The restframe equivalent width 
of Si IIA1808, Wr = 1.13 A is also the largest ever de- 
tected in any extragalactic sightline. Assuming the latter 
line is located on the linear part of the curve of growth 
implies that logA^(Si II) = 16.3. This is a very con- 
servative lower limit to the actual column density and 
translates into a conservative lower limit on the metallic- 
ity of [Si/H] > -0.8. Note also that Si could be depleted 
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TABLE 3 

Ionic column densities of the GRB system. 



Ion 


Transition lines used 


log TV 




(A) 


(cm~^) 


Si II 


1808 


16.6(0)±0.3{2) 


Znll 


2026; 2062 


13.6(7)±0.6{5) 


Mgl 


2026 


14.4(1)±0.6{4) 


CI 


1560; 1656 


14.9(5)±0.2{5) 


Fell 


2249; 2260=" 


16.1(6)±0.1{8) 


Cr II 


2056; 2060; 2066=^ 


14.0(4)±0.4{0) 


Ni II 


1741 


14.8(9)±0.2(8) 


Mn II 


2576 


13.8(9)±0.1(8) 



Note. — The broadening parameter was b =123±74 km s ^ in 
all cases. 

onto dust grains (jPetitiean et al.ll2002f ). 

We further performed Voigt-profile fitting (see Table[31) 
using some of the weakest (least saturated) metal lines 
(see Fig. |3]). The absorption hne at Aobs ~ 7000 A 
is likely to be a blend of Zn II and Mg IA2026 with 
possibly similar strengths. It is therefore important to 
fit both transitions simultaneously. The other Zn II 
(A2062) line is also blended (with one of the Cr II triplet 
lines). Therefore, the best line to estimate the metal- 
licity remains Si IIA1808. All lines were simultaneously 
fitted with a single component constraining their red- 
shifts as well as their broadening parameter (assumed 
to be purely turbulent) to have identical values. We 
find [Si/H] = —0.46 ± 0.38 and also consistently get 
[Zn/H] = -0.50 ± 0.68. Assuming the Zn IIA2026 line 
is located on the linear part of the curve-of-growth im- 
plies that [Zn/H]i-0.35 but the actual metallicity could 
be lower than that due to the unknown contribution of 
the Mg IA2026 line to the observed feature (see above). 

The large error bars on the measured column densities 
are mainly a consequence of the large uncertainty on the 
broadening parameter b which we find to be: b = I23±74 
km s~^. One might argue that b could be significantly 
smaller (e.g., 6 < 30 km s^^) in which case the col- 
umn densities would be much larger than found in our 
best fit. However, it has been shown that the width of 
metal line profiles in DLAs is larger at higher metallicities 
(|Ledoux et al.ll2006HProchaska et alT 20081. Therefore, a 
large b value is consistent with the relatively high metal- 
licity infer red, which is amon g the highest measured for 
any GRB (|Fvnbo et al. l [200l . 

The measured abundance ratios are consistent with 
what is seen in other GRB- or QSO-DLAs. We find that 
the mean depletion of Fe-peak elements (Fe, Cr and Mn) 
compared to Si is ~ 0.5 dex. This is not particularly 
high, but yet significant implying that dust is present in 
this system. According to our best fit, the Si metallicity 
could be as large as solar. This in fact is quite likely 
taking into account that Si is affected by dust depletion 
fi'etitiean et al . 2002). 

2.1.2. Neutral species 

In the afterglow spectrum, we clearly identify neutral 
carbon lines, namely C I AAI560,1656 (blended with C I* 
and C I**; see Fig. [3]). These are the first ever clear detec- 
tions of C I in a GRB-DLA which typically have very low 
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Fig. 3. — Our Voigt-profile fits to low-ionization metal lines. 
Best-fitting parameter values are given in Table |3] 

C I/C II ratios (jProchaska et al.l[2007l ). The two features 
are very strong. We measured Wr = 1.4 A at Aobs ~ 5720 
A which is dominated by C I A1656 absorption. Using the 
above fit results, we get log N(C I)/A^(Si II) = -1.65. 
This can be conipare d to QSO-DLAs from Fig. 10 of 
ISrianand et al.l ()2005D showing that II2 is most likely 
present as well in this system, although this can not be 
verified as it falls too much in the blue into the Lya 
forest, and because of low signal to noise and spectral 
resolution. On the other hand, the strength of the de- 
tected C IV doublet line is only modest and similar to 
most other GRB sight-lines. 

2.1.3. Non detection of Fe If m GRB 070802. 

Over the past few years, excited lines of Fe II 
have been detected along several GR B sightlines (e.g., 
iSavadio fc Faia[200l[Chen et al.ll2005l) . which have been 
shown to be due to indirect exci tation by UV photons 
p roduced by the GRB afterglow fProchaska et al."2006!; 
Dessauges-Zavadsky ct al. 2006; Vrccswiik et al. 2007). 
For two specific GRBs, modeling of the afterglow flux 
exciting a cloud of Fe II and/or Ni II atoms has led to 
an estimate of the distance between the GRB and the 
bulk of the neutral absorption material (responsible for 
low-ionization lines from ions such as C I, O I, Fe II, 
Cr II, Mn II, Si II, Zn II, Ni II): 2 kpc for GRB 060418 
(Vreeswijk et al. 2007), and 0.5 kpc for GRB 050730 
(Ledoux et al. 2009, submitted). For the GRB 070802 
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sightline we detect consistent features at the expected 
positions of the Fe II*A2396 and Fe II**A2405 Unes at 
z=2.455 but, since the hues are hkely to be at the same 
time shghtly saturated and blended, we cannot infer from 
them secure column densities. Based on the distance es- 
timates in previous GRBs we note that the bulk of the 
neutral material, which likely corresponds to the region 
where the extinction bump is originating, is probably at 
least 0.5 kpc away from GRB 070802. 

2.1.4. Constraints on the z — 2.078 foreground absorber 
extinction 

The foreground absorption system at z = 2.078 ex- 
hibits strong Mg II (W^^^yge ^ 3 9 ^-j g^^^^ II lines 
with equivalent width ratio W^^^2796/^A2600 < 2. Ad- 
ditionally the detected Mg IA2852 line has Wr = 1.0 
A. These charac teristics make th is absorber very likely a 
DLA (see, e.g., iRao et al.|[200l their fig. 11). We note 
that such a strong Mg II absorber is very rare for QSO 
sightlines, but it has been suggested that GRB sightlines 
may have strong er absorbers although the reason for this 
is still unknown (jProchter et al.|[2006[ ) . 

To constrain Ay from the foreground absorber we look 
at the metal lines as metallicities and dust depiction 
factors (and h ence the dust content) correlate in DLAs 
(jLedoux et al.| [2003). The non-detection of Zn II A2026 
leads to an upper limit on the zinc column density. Un- 
fortunately, the constraint we can derive from the spec- 
trum, logA^(Zn II) < 13.5 at the 3tT level, is fairly weak 
because there is an absorption feature immediately blue- 
wards of the expected position of the Zn II line. More- 
over, we cannot use Si II instead of Zn II because the 
expected position of the Si IIA1808 line at z = 2.078 is 
blended with Fe II A 1608 from the GRB host galaxy ab- 
sorber aX z — 2.455. The best constraint on Ay therefore 
comes from the observed Fe II lines at z = 2.078. The 
ratio of the weakest detected Fe II line, i.e., Fe IIA2374, 
at z = 2.455 to that at z = 2.078 is 4.3. Because 
of saturation effects, this value is a lower limit on the 
ratio of the Fe II column densities. Altogether, the 
weaker Fe II lines at z = 2.078 and the non-detection 
of Zn IIA2026 at z = 2.078 are an indication of lower 
metallicity because in DLAs the dust depletion factor 
is only a slowly decreasing function of met allicity (e.g., 
iLedoux et alll2003l : iNoterdaeme et al.ll2008f) . A factor of 
> 4.3 leads tologjy(Zn I I) < 13 and therefore Ay < 0.25 
(iVladilo & Perouxl[2005L their fig. 1). 

An alternative way to estimate the contribution of the 
foreground absorbe r is to us the relation reported by 
iMenard et all ()2008l . their eg. (18)) which uses the equiv- 
alent width of the Mg II line to estimate E{B — V). As- 
suming a standard MW extinction law, with Ry = 3.1, 
this translates into A{V) w 0.06 which is well below the 
limit derived above. 

2.2. X-ray observations 

The X-ray telescope onboard Swift observed the af- 
terglow of GRB 070802 beginning observations in pho- 
ton counting mode at 147 s after the trigger time {to = 
07:07:26 on 2 August 2007). The data were extracted and 
reduced in a standard way using the HEAsoft software 
(version 6.2) and the most recent calibration files. The X- 
ray spectrum (with a total exposure time of 3419 s in the 
interval to + 147-^0 + 16800 s) was extracted and fit with 




Fig. 4.— The field of GRB 070802. The left panel shows the 
afterglow as detected about 1 hr after the burst. The right panel 
shows a deep late-time exposure of the same field (65 days after the 
burst), taken when the afterglow had faded away, and revealing a 
host galaxy at the GRB position. Each panel is 30" X 30" wide. 

a power-law with absorption fixed at the G alactic level 
(2.9 X 102°cm-2) (jDickev fc LockmanI [T990l) and freely 
variable absorption at the redshift of the host galaxy. 
The best fit resulted in a power-law with a photon spec- 
tral index r = 2.02+^:^^ (/3 = 1.02+^;^^, 68% confidence 
level) and an equivalent absorbing hydrogen column den- 
sity of 9.7^0 4 X 10^^ cm~^ at z = 2.45 assuming solar 
abundances. To produce the X-ray segment of the near- 
infrared to X-ray spectral energy distribution (SED), this 
X-ray spectrum was normalized to the fiux level obtained 
from the X-ray lightcurve at the SED time {to + 3517 s) 
and corrected for both absorptions. No significant evi- 
dence (above the 1 a level) was found for spectral varia- 
tions over the period of the spectrum, by comparing the 
first 400 s and the remainder of the observation. 

2.3. The host galaxy 

K- and i?-band observations of the host galaxy of 
GRB 070802 were obtained with the ESQ VLT during 
the nights of 27 September 2007 and 5 October 2007 
(57 and 65 days after the GRB). At the afterglow posi- 
tion, a source is clearly detected with R = 25.03 ± 0.10 
(Fig.|4]). A faint source is also visible in K, with magni- 
tude K = 21.70 ±0.25. The centroid of the host is offset 
by 0.15 ± 0.04 arcsec with respect to the afterglow posi- 
tion (1.2 kpc at z = 2.455). From the observed K mag- 
nitude (rest frame 6270 A), and assuming a power-law 
spectrum consistent with the R — K color {f^, oc 
we can estimate the absolute i3-band magnitude of the 
host to be Mb — —21.0 (uncorrected for extinction) 
which is fairly luminous f or a GRB host (jFruchter et al.l 
l200llSavaglio et al.ll2009l) . 

3. THE EXTINCTION CURVE 

The spectrum shows strong extinction features, with a 
clear detection of a bump near A = 2175 A in the rest- 
frame of the host (Fig.[T|). To extract the continuum part 
of the spectrum to use for the extinction curve fitting we 
start by selecting the wavelength range 4500-9000A. We 
then remove emission and absorption lines by recursively 
fitting the spectrum to an 8th order polynomial removing 
parts which deviate by more than 2.2a (see Fig. [1]). 

The intrinsic fiux of the afterglow is modeled as 

rr = f.o{^y' (1) 

(where v is frequency, (3 the intrinsic slope and f^^o is the 
flux at i^o) but is partially extinguished due to dust along 
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the line of sight. The observed flux is therefore given by 

^obs^^i„t.iO-0.4Ar-- (2) 

= /.,o(^)''lO-"-<^^+^?^') (3) 

where the total extinction is A^"**"' = Ax + Af'^^ Ax is 
the extragalactic extinction along the line of sight as a 
function of the wavelength A and A^^' is the extinction 
in the Milky Way. The extragalactic extinction can be 
caused by both extinction in the host of the GRB and in 
foreground objects, i.e. Ax = A\°^^+A^^''°. The observed 
flux is corrected for extinction in the Milky Way using 
the maps of Schlcgcl ct al. (.1998|, E{B - V) = 0.025, 
Rv = 3.1), leaving 

/. = /.,o(^) \o-'-''^^ (4) 

which is the intrinsic flux of the burst affected only by 
extragalactic extinction along the line of sight. Solving 
for the extinction, we find 

where x = 1/A = v/c is the wavenumber. The abso- 
lute extinction along the line of sight can therefore be 
determined if the underlying spectral slope f3 and the 
normalization f^ Q at xq are known. 

The simple shape of the intrinsic spectra of GRB af- 
terglows (see § II. 3p allows us to constrain both /3 and the 
normalization. Assuming that there is no cooling break 
in the intrinsic spectrum, one may use the slope and 
normalization derived from the X-ray data, i.e. (3 = (3x- 
Alternatively, assuming that there is a cooling break in 
the spectrum, the slope must be /3 = /3x — 0.5 while 
the normalization depends on the location of the cool- 
ing break. These two possibilities will both be addressed 
here and are depicted in Fig.[5l In addition, we will start 
by assuming that all the extinction is caused by dust in 
the host galaxy of GRB 070802, i.e. that 

Ax=A\°'\ (7) 

The possible contribution from the foreground Mg II ab- 
sorption systems (§ I2.1.4p will be addressed in § 13.31 

We fit the derived extinction curves to four types of ex- 
tin ction laws, i.e. MW type extinction as parametrized 
bv ICardelli et all (|1989f ). LMC and S MC type extinc- 
tion as parametrized by iPeil (1X9921) and the modi- 
fied p arametrization for the UV of iFitzpatrick fc Massal 
()2007l FM). A d escrip tion of these extinction laws 
and how the | Pei| (Il992f) parametrization compares to 
the I Gordon et all ()2003D SMC and LMC extinction 
cu rves is given in Appendix [Al The parametrization 
of IFitzpatrick fc Massal (|2007| ) does not directly tie the 
slope to the bump, but keeps the parameters of each 
independent. It therefore has much greater freedom in 
fitting the curve and is able to accurately trace the de- 
rived extinction. The parametrization however assumes 
that Ay and Ry have been independently derived from 
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Fig. 5. — Spectral energy distribution of the afterglow of 
GRB 070802 at to + 3517 s. The plot shows X-ray data (trian- 
gles, corrected to the time frame of the optical photometric points) 
and the derived intrinsic powerlaw spectrum (black solid line) and 
the 68% confidence levels (black dashed lines). The spectrum (red 
curve) is scaled to the i?-band. The squares are the photometric 
points from this pap er, the diamonds are the photometric points of 
IKriihler et al.l I I2008I ). Assuming that the intrinsic slope of the GRB 
can be described as a single powerlaw (i.e. (3 = Px = l-02^Q'|[g), 
the difference in the observed spectrum and the extrapolated X-ray 
spectrum is interpreted as being due to extinction. The blue lines 
correspond to the upper and lower limits of the intrinsic power law 
given a cooling break in the spectrum. The intrinsic slope in the 
optical is then /3 = f3x — 0.5 = 0.52^Q'jg and the normalization 
is determined by the location of the cooling break, which can be 
anywhere between the optical and X-ray data sets with the addi- 
tional constraint that the optical data points can not be brighter 
than the intrinsic curve. The intrinsic powerlaw spectrum with a 
cooling break corresponding to the 68% confidence levels have been 
omitted for clarity. 

the optical extinction. As our data do not reach into the 
restframe optical, we have opted to fix Ry = 3.1 and 
allow Ay to vary while fitting the UV. We also tested 
choosing Ry = 2.7 and Ry = 3.5 and find that it leads 
to a change in Ay and ci (one of the parameters of the 
FM extinction law, see App. |^ of around 10% but does 
not significantly affect the other parameters. 

We compare the ext rapolation of our fits to the infrared 
photometric points of lKriihler et al.1 (|2008D . We used the 
spectral energy distribution (SED) from GROND cre- 
ated by citetkruhlerOS to derive the near-infrared data 
we sho w here simply by ex trapolating the mean SED 
time of lKriihler et all (|2008f ) to our mean SED time us- 
ing our estimate of the optical decay rate. We note that 
these points are not used in constraining the fits. 

3.1. The extinction assuming no cooling break 

If there is no cooling break between the X-ray and opti- 
cal regimes, clearly the X-ray and optical emission lie on 
the same power-law. Under this assumption, the intrin- 
sic slope and normalisation in the optical can therefore 
be derived from a simple extrapolation of the X-ray data. 
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Any discrepancy between the observed optical flux and 
the extrapolation is then interpreted as extinction along 
the line of sight. A slight complication is introduced 
by the non-simultaneity of the observations in different 
bands, since the flux decays with time. However the un- 
certainty introduced by this correction is not large, since 
the optical photometry mean times are I — to + 3517 s, 
V = to + 3627 s and R is interpolated between observa- 
tions at to + 3249 s and to — 3949 s, and the X-ray flux is 
derived from a lightcurve that covers this time period. 

Fig. [5] shows the X-ray spectrum (scaled to the SED 
time, to + 3517) and its extrapolation into the optical. 
The fluxed optical spectrum (mean time tg + 2.2 hours) 
was scaled to the i?-band data point. It is worth not- 
ing that the scaled optical spectrum agrees with the ob- 
served V and / photometric data, indicating that the 
small temporal extrapolations are correct and that there 
is no significant spectral change in the optical in this pe- 
riod. The extinction curve, shown in the top panel of 
Fig. [6l is then derived from the difference between the 
spectrum and the extrapolated slope. 

To study the properties of the extinction curve, we 
fit it to the four previously mentioned extinction laws 
with the additional constraint that the extinction must 
go to zero for infinite wavelengths. The derived extinc- 
tion and the fits are shown in the top panel of Fig. [6] 
and the parameters of the fit are given in Tables |4] and 
[51 The freedom of the FM parametrization allows for a 
good fit (xVd-O.f. = 261/968), the LMC also provides a 
reasonable fit (x^/d.o.f. = 815/975), while the SMC and 
MW fail to reproduce the extinct ion curve. However , 
the infrared photometry points of iKriihler et al.l ()2008f ) 
clearly do not agree with any of the fits and require the 
extinction curve to take a very sharp break in its steep- 
ness at X < 1.5 /im~^ which would suggest that the host 
of GRB 070802 has a very strange form of extinction 
law. However, a more plausible explanation is that the 
assumed intrinsic slope is wrong, either due to the ex- 
trapolation or because of the presence of a cooling break 
in the spectrum (see ij 13. 2p . 

From Fig. [5] we see that due to the several order of 
magnitudes extrapolation, the 1 a difference in f3 trans- 
lates into a difference of ~1 magnitude in the optical. 
While going towards steeper slopes will make the appar- 
ent 'break' in the extinction curve greater, a shallower 
slope will act to reduce it as seen in Fig. [6l Redoing 
the analysis for a shallower slope corresponding to 1 cr 
deviation (i.e. (3 = 0.87) results in the extinction curve 
shown in the middle panel of Fig. [6l The parameters 
of the fit, for both 1 a deviations are given in Tables |4] 
and El The goodness of fit for the MW, SMC and LMC 
is improved for the shallower slope, while the effect on 
the FM fit is minimal. The fits are also more consistent 
with the addit i onal p hotometric points in the infrared of 
IKriihler et all (|2008f ). We note however th at most GRBs 
have /3x > 1 (see e.g.. iLiang et al.l [20071 their fig. 3, 
with T = P+l). 




Fig. 6. — The derived extinction curve (black line) of the 
GRB 070802 host galaxy. The curve has been smoothed for clarity 
in presentation. The curve has been fit by an SMC (blue long- 
dashed line), LMC (cyan dotted line), MW with Ry fixed (red 
solid line) and allowed to vary (red dashed line) and finally using 
the FM parametrization (yellow dash-dot line). The filled dia- 
monds are the photometric points presented in this paper used to 
scale the spectrum while the empty diamonds are the photometric 
points of Kriihlcr et al. (2008). Note that the fits were done using 
only the spectroscopic data. Top panel: The derived extinction 
curve assuming no cooling break between the optical and X-rays 
in the intrinsic spectrum of the burst, i.e. 13 = fix = 1.02. The 
best fit is obtained using the FM parametrization. Of the three lo- 
cal extinction laws (MW, SMC, LMC), the LMC provides the best 
fit. The photometric points in the infrared are in poor agreement 
with the extrapolations of the fits and would require a break in 
the extinction curve at x < 1.5 /im~^. Middle panel: The same 
as the top panel, but using the 1 tr deviation value for /3 towards 
a shallower intrinsic slope, /3 = 0.87. In this case, the photometric 
points are in better agreement with the predicted extinction curve 
and the sharp break in steepness at around a; < 1.5 /im~^ is not 
required. However, it implies fix < 1- Bottom panel: The de- 
rived extinction curve (black line) assuming a cooling break in the 
intrinsic spectrum, i.e. (3 = Px — 0.5 = 0.52, and shifted by a con- 
stant Aac which depends on the location of the cooling break. The 
FM parametrization gives the best fit to the data while the LMC 
provides the best fit a moung the thr e e loca l extinction laws. The 
photometric points of Kri ihler et al.l II2008I ) are inconsistent with 
the extrapolation of the SMC and MW (with a fixed Ry = 3.1), 
but are in agreement with the other fits. 



3.2. The extinction curve assuming a cooling break 

A cooling break in the intrinsic spectrum occurring be- 
tween the X-ray and the optical data, i.e. /3 — (3x — 0.5, 
would result in an intrinsic spectrum at optical wave- 
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TABLE 4 

Parameters of the fits with no cooling break 



/3 


Type 






y1 fore 
V 


/dof 


1.02 


MW 


0.87+0.03 


1.8+0.1 




2281/974 


1.02 


MW 


1.533+0.002 


3.1 




2755/975 


1.02 


FM 


1.341+0.002 


3.1 




261/968 


1.02 


LMC 


1.474+0.002 






815/975 


1.02 


SMC 


1.269+0.002 






7848/975 


1.02 


MW+SMC 


1.274 + 0.002 


3.1 


(0.25) 


1102/974 


1.02 


LMC+SMC 


1.474 + 0.002 




(0) 


815/974 


1.19 


MW 


1.92+0.03 


2.9+0.1 




4670/974 


1.19 


MW 


2.052+0.002 


3.1 




4690/975 


1.19 


FM 


1.820+0.002 


3.1 




321/968 


1.19 


LMC 


1.976+0.002 






2654/975 


1.19 


SMC 


1.705+0.002 






16698/975 


1.19 


MW+SMC 


1.795 + 0.002 


3.1 


(0.25) 


2821/974 


1.19 


LMC+SMC 


1.976 + 0.002 




(0) 


2654/974 


0.87 


MW 


0.32+0.03 


0.9+0.1 




1229/974 


0.87 


MW 


1.052+0.002 


3.1 




1859/975 


0.87 


FM 


0.905+0.002 


3.1 




287/968 


0.87 


LMC 


1.018+0.002 






338/975 


0.87 


SMC 


0.883+0.002 






2692/975 


0.87 


MW+SMC 


0.796 + 0.002 


3.1 


(0.25) 


649/974 


0.87 


LMC+SMC 


1.018 + 0.002 




(0) 


338/974 



Note. — The parameters of the fits for the different extinction 
laws assuming no cooling break between the optical and X-rays in 
the spectrum. The slope in the optical is then jS = fix = l-02^g :[g. 
Numbers in italics were kept fixed, while numbers in brackets 
reached the limits of their allowed range. The error bars are the 
formal 1 cr errors from the minimization. For the fits assuming 



a foreground contribution (at z = 2.08), Ay is the extinction of 
the host and A^''° is the extinction of the foreground object, both 
in their respective restframes. 



lengths of the form 



j'intr r 

Ju — Jv.» 



(8) 



where the normalization f^ o depends on the location of 
the cooling break, Xbroak, by 



■^break — -^0 



The absolute extinction is then given by 





-2.5 log. 



-2.51ogio - 



(9) 

(10) 
(11) 



The first term of the sum can be determined as before 
from the X-ray data, while the second part of the sum 

^ 2.5 log,, ^ = 2.51og,n (^)° ' (12) 



=10 7^ = 2.51ogio 



is an additional unknown to be added to the fits. Its 
effect is to shift the extinction curve up or down (i.e., 
the value of A\ + A^o at a: = is Aoo)- Constrain- 
ing the break to occur between the optical and X-ray 
regimes translates into constraints on the allowed values 
of Aoo corresponding to setting Xbrcak to be the upper 
X of the optical data and the lower x of the X-ray data. 




Fig. 7. — The effect on an "observed" extinction curve assuming 
that the extinction is caused by a MW type extinction in addition 
to an SMC foreground object (located a.t z = 2.45 and z = 2.08 
respectively). The three lines correspond to Ay = 1.0, Ay'^° = 0.0 
(solid line). Ay = 0.75, A^"''" = 0.25 (dotted line) and Ay = 0.5, 
Ay'^° = 0.5 (dashed line) where the Ay values are quoted in their 
restframes. 

An additional requirement is that the extinction remains 
positive. 

We fit the three local extinction laws (MW, SMC and 
LMC) to the curve in addition to the FM parametriza- 
tion. Again we find that of the three local extinction laws 
the LMC provides the best fit (xVd.o.f. = 311/974), 
but the FM parametrization still gives the best overall 
fit (x^/d.o.f. = 253/96 7). In this ins t ance, the infrared 
photometric points of IKrtihler et al.l ()2008[ ) agree with 
the extrapolated LMC and FM fits. 

3.3. Possible contribution from a foreground object 

The MW extinction law provides a poor fit to the ob- 
served extinction, in part because the observed bump is 
shallower than the predicted MW bump and the rise into 
the UV is steeper (see e.g. Fig. [5]). However, as Fig. [7] 
shows, such effects can be caused if part of the extinction 
is due to a foreground object without a bump. 

The spectrum of GRB 070802 shows two intervening 
Mg II absorption systems which might also cause part 
of the reddening. Neither foreground system shows any 
signs of a 2175 A bump. If the foreground absorber had 
a 2175 A bump in its restframe, it would act to broaden 
and shift the peak of the observed bump in the restframe 
of the host, which is inconsistent with the data. There- 
fore, the most straightforward way to model them is with 
an SMC type extinction law. As this extinction law is 
roughly linear, there is no way of separating the effect of 
two systems without additional constraints on the dust 
content of the two systems. As the absorber at z — 2.29 
has much weaker absorption lines, and therefore by as- 
sumption much weaker extinction, we will only look at 
the possible contribution from the system at z = 2.078. 

We will therefore limit ourselves to studying the pos- 
sibility of the extinction being of the form 



Ax^A 



host 



fore 



(13) 



where A\°^^ 



is the extinction in the host galaxy and A^^'^° 
is the extinction in the foreground Mg II absorber. We 
have already seen that the extinction curve is poorly fit 
by an SMC type extinction, as it lacks the characteristic 
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TABLE 5 

Additional parameters of the FM fits 



/3 


Cl 


C2 


C3 


C4 


C5 




7 


1.02 


0.080 ±0.001 


1.025 ±0.001 


2.747 ± 0.007 


0.355 ± 0.003 


5.239 ±0.013 


4.583 ± 0.004 


1.084 ± 0.002 


1.19 


0.072 ± 0.001 


1.033 ±0.001 


2.719 ±0.005 


0.278 ± 0.005 


5.889 ±0.031 


4.463 ± 0.003 


1.134 ±0.002 


0.87 


0.079 ± 0.001 


1.011 ±0.002 


2.823 ±0.010 


0.531 ± 0.005 


4.779 ±0.012 


4.617 ±0.007 


1.040 ± 0.003 


0.52 


0.083 ± 0.001 


1.018 ±0.002 


2.683 ±0.007 


0.317 ±0.002 


4.672 ±0.010 


4.642 ± 0.005 


1.094 ± 0.002 



0.69 
0.37 



0.071 ± 0.003 
0.073 ± 0.001 



1.027 ±0.001 
1.023 ±0.002 



2.723 ±0.008 
2.837 ±0.012 



0.366 ± 0.002 
0.440 ± 0.003 



4.847 ±0.012 
4.160 ±0.011 



4.606 ± 0.005 
4.667 ± 0.008 



1.095 ±0.002 
1.065 ± 0.004 



6 
5 

B 4 

CO 

£ 3 
< 2 
1 



MW+SMC 

VIC+SMC 





MW+SMC 

LMC+SMC 



Note. — The additional parameters of the FM parametrization for fits with or without a cooling break. The errors quoted are the formal 
Iff errors. The uncertainty in the parameters is dominated by the uncertainty in /3. The top three lines correspond the /3 = /3x = 1.02 and 
the 1 ff deviations, while the lower three lines correspond to /3 = I3x — 0.5 = 0.52 and the 1 a deviations. 

bump. We will therefore parametrize the extinction of 
the host galaxy of GRB 070802 as either a MW (with 
fixed Rv — 3.1) or an LMC type, while the foreground 
absorber will be parametrized as an SMC type. In addi- 
tion, we place the limit Ay'^° < 0.25 for the foreground 
absorber based on the discussion in § 12.1.41 

The results are shown in Fig. [5] and tabulated in Ta- 
bles Hand H We find that for the LMC, the addition of 
the foreground contributor does not improve the fits and 
the best fits are found by setting Ay'^'^ to be zero. For 
the MW, the fits are improved and the best fit is found 
by setting Ay'^° to its maximum value of 0.25. The MW 
fits are still worse than the LMC leading us to conclude 
that the extinction of GRB 070802 is not well fitted by 
a MW type extinction. 

3.4. Summary of the extinction curve properties 

Fig. [9] shows the scaled extinction curve for 
GRB 070802 for /3 = 0.87. We clearly detect the 2175 A 
bump in the extinction curve of GRB 070802. The data 
sample the bump very well and its detection does not de- 
pend on whether we assume cooling breaks or not in the 
spectrum of the GRB. It is one of the most robust detec- 
tions of the bump in extragalactic environments to date 
and currently the highest redshift detection at z=2.45 
(corresponding to 2.5 Gyr assuming a flat universe with 
flm = 0.3, Qa = 0.7 and Hq = 73 km s'^ Mpc'^). 

For the three local extinction laws, the LMC provides 
the best fit (i.e., it has the lowest per degree of freedom 
(d.o.f.)) to the shape of the derived extinction curve, re- 
gardless of whether we assume there is a cooling break 
or not between the X-rays and the optical, and whether 
we take into account a possible contribution from the 
foreground absorbers or not (see Tables 2] and [S]) . The 
SMC clearly provides a poor fit in all cases as the bump 
is completely missing, while the bump of the MW ex- 
tinction law is too strong and the rise into the UV not 
steep enough. A foreground absorber with SMC type 
extinction could make the bump shallower and the UV 
rise steeper (see Fig. [7]), however, the absorption would 
need to be larger than the upper limits placed on Ay"'^ 
in §[2X1 

The best fits are obtained using the FM parametriza- 
tion. This is not surprising, as it has more freedom in 
tracing the shape of the curve (see Appendix [X| . The 
three parameters giving the shape of the bump arc only 
weakly dependent on the assumed (3 and whether we as- 
sume a cooling break or not (Table . The bump is 
found to be centered at Xc ~ 4.6 ± 0.1 /^m~^, the width 




2 4 6 8 

Fig. 8. — Extinction curve fits taking into account a possible 
contribution from the foreground absorber at z = 2.078. The ex- 
tinction of the host is taken to be a MW (red solid line) or an LMC 
(cyan dotted line). The extinction of the foreground absorber is 
taken to be an SMC (as there is no evidence of a bump in the data) 
with an upper Umit of A^"''" < 0.25 (see q |2X4t . The filled dia- 
monds are the photometric points presented in this paper used to 
scale the spectrum while the empty diamonds are the photometric 
points of Kriihler ct al. (2008). Note that the fits were done us- 
ing only the spectroscopic data. The top panel shows the derived 
extinction assuming no cooling break. The middle panel shows 
the 1 CT deviation towards a shallower intrinsic slope. The bottom 
panel shows the derived extinction assuming a cooling break. For 
the LMC fits the contribution of the foreground extinction is not 
found to improve the fits (setting Ay'^^ = 0). The MW fits are 
improved by adding an SMC contribution, however, they are still 
worse than LMC only, and would require Ay'^'^ > 0.25. 
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TABLE 6 

Parameters of the fits assuming a cooling break 






TvDe 


Av 




Rv 


A tore 


/dof 


0.52 


MW 


(Q) 


— 1.47+0.04 








0.52 


MW 


0.68+0.02 


— 1.14+0.06 


3.1 




2517/974 


0.52 


FM 


1 2,59+0 002 


f-S 2~) 


3.1 




253/967 


0.52 


LMC 


1.27+0.02 


-2.86+0.07 






311/974 


0.52 


SMC 


0.52+0.01 


-0.97+0.04 






1078/974 


0.52 


MW+SMC 


0.69 + 0.02 


-1.85 + 0.06 


3.1 


(0.25) 


1030/973 


0.52 


LMC+SMC 


1.27 + 0.02 


-2.86 + 0.07 




(0) 


311/973 


0.69 


MW 




(-0 5~) 


0+0.01 




994/973 

\J 1 <J i '-J 


0.69 


MW 


964+0 002 


(-0 5~) 


3.1 




2155/974 


0.69 


FM 


1.25+0.03 


— 1.76+0.08 


3.1 




253/967 


0. 69 


LMC 


1.17+0.02 


-1.17+0.06 






307/974 


0.69 


SMC 


0.807+0.002 


(-0.5) 






1956/974 


0.69 


MW+SMC 


0.706 + 0.002 


(-0.5) 


3.1 


(0.25) 


784/973 


0.69 


LMC+SMC 


1.17 + 0.02 


-1.17 + 0.06 




(0) 


307/973 


0.37 


MW 


(0) 


-2.67+0.04 


+ 0.005 




1125/973 


0.37 


MW 


0.73+0.02 


-2.57+0.06 


3.1 




2671/974 


0.37 


FM 


0.807+0.002 


(-3.2) 


3.1 




410/967 


0.37 


LMC 


0.937+0.002 


(-3.2) 






559/974 


0.37 


SMC 


0.58+0.01 


-2.45+0.04 






1060/974 


0.37 


MW+SMC 


0.69 + 0.02 


-3.15 + 0.06 


3.1 


(0.25) 


1184/973 


0.37 


LMC+SMC 


0.78 + 0.02 


(-3.2) 




0.16 + 0.02 


512/973 



Note. — The parameters of the fits for the different extinction laws assuming a cooling break between the optical and the X-rays in the 
spectrum. The slope in the optical is then /3 = fix — 0.5 = 0.52^Q'j^g. Numbers in italics were kept fixed, while numbers in brackets reached 
the limits of their allowed range. The error bars are the formal 1 a errors from the ^ minimization. For the fits assuming a foreground 
contribution (at z 
restframes. 



2.08), Ay is the extinction of the host and A^y^'^ is the extinction of the foreground object, both in their respecti- 



of the bump is found to be 7 « 1.08 ± 0.05 [lvoT^ while 
the 'strength' (i.e. its height above the linear extinc- 
tion, see Appendix]^ is C3 w 2.7 ± 0.1. Therefore, both 
the area of the bump Abump = ''i'C3/(27) and its maxi- 
mum height above the linear extinction i?bump = c^/l"^ 
are well defined. The relative strength of the bump 
^bum p I Ay = C3/(7^i?y) (as defined by iGordon et al.l 
I2OO3I) is also well defined, although here the uncertainty 
is dominated by the uncertainty in Ay (or equivalently 
Rv)- We note that the valu e of C3 is the s ame a s the 
average C3 = 2.7 ± 0.1 that iGordon et all (|2003D find 
for the LMC average sample while the width 7 is a bit 
wider compared to their 7 = 0.93 ±0.02 ^m~^ (although 
it is still within their scatter). The value of C3 is how- 
ever higher than IGordon et all (|2003f ) find for the LMC2 
sample (03 = 1.5 ± 0.1) although it is within the scatter. 

We estimate the amount of dust extinction to be 
Av = 1.341 ± 0.002 for /3 = 1.02 and Av = 1.259 ± 0.002 
for (3 = 0.52 given the FM parametrization. The LMC 
fits correspondingly give Av = 1.474±0.002 for /3 = 1.02 
and Av = 1.27 ± 0.02 for (3 = 0.52. Taking into account 
the possible 1 a deviation of the slope (excluding steeper 
slopes than fi = 1.02, as this would be in strong dis- 
agree ment with the photometric points of iKriihler et al.l 
12008). we estimate Av — 0.8-1.5 in the host along the 
fine of sight to GRB 070802. 

4. DISCUSSION 

The nature of the interstellar extinction peak at 2175 A 
remains poorly understood more than 40 years after 
its discovery by Stecher (1965). As its detection has 
been limited to the Milky Way with only a few excep- 
tions it has proven hard to search for correlations be- 
tween the dust environment and the detection or non- 



detection of the bump. The detection of the bump in 
the spectrum of GRB 070802 is interesting in itself for 
two reasons, i.e., as being the highest redshift detection 
of the bump and as being the first robust detection of 
the bump in a GRB host galaxy. It shows that the 
carrier of the 2175 A bump, which is characteristic for 
Milky Way type dust, was in place 2.6 Gyr after the Big 
Bang (when the Universe was only 20 % of its current 
age). It also shows that the conditions for both forming 
and not destroying the 2175 A were satisfied in a GRB 
host galaxy - surprising in view of the fact that most 
GRB hos t galaxies are faint, blue, young, low-metallicit: 
galaxi es (iLe Floc'h et all l2003t IChristensen et al.l 1200 



Fruch ter et al.l 



2006ir in contrast to the massive, evolved 



and chemically enriched Milky Way. 

Moreover, the detection of 2 175 A as well as a very 
rich spectrum of redshifted UV metal absorption lines 
allows us to explore various hypotheses for the origin of 
the bump - an experiment that cannot easily be done 
along lines of sight in the Milky Way. Below we ex- 
plore whether there are any other properties of the galaxy 
which are correlated with the presence of the bump. We 
discuss what separates GRB 070802 from other GRBs 
with reddening but no detected bump and compare the 
properties of the host galaxy with other galaxies for 
which the extinction curve has been determined. 



4.1. Metallicity 

As mentioned in § II. 2[ based on evidence from the 
MW, LMC and SMC, the simplest hypothesis would 
be that the strength of the 2175 A b ump is simply 
contro lled by the metallicity. However, ISavaglio et alj 
(|200l infer a metalhcity of [Zn/H]= -0.13 ± 0.25 for 
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2 4 6 8 

1/X (|inT^) 



Fig. 9. — An absolute extinction curve for the afterglow GRB 070802 at 2 = 2.45. It is based on the simplest model of the afterglow 
that is consistent with all the available data (single power-law spectrum, /3 = 0.87, corresponding to a 1 cr deviation from the best fit) and 
with the constraint that the extinction should be zero at 1/A = 0. The extinction curve shows a clear bump at 2175 A. Also shown are the 
best fits (to the spectroscopic data) for an LMC extinction (dashed red line) and MW extinction (dotted blue line), with the LMC clearly 
providing the better fit. To obtain the absolute extinction, Av was derived using a linear interpolation of the two datapoints closest to 
the restframe y-band. The filled diamo nds are the ph otometric points presented in this paper used to scale the spectrum while the empty 
diamonds are the photometric points oflK riihler et al. (2008). 



GRB 000926 (compared to [Zn/H]= -0.50 ± 0.68 we de- 
rive for GRB 070802) and the H I column density for 
GRB 000926 is similar to that of GRB 070802. This 
would indicate that the extinction curve of GRB 000926 
should also have a bump if metallicity was its only driver. 
However, the extinction curve derived for GRB 000926 is 
inconsistent with those of the MW and LMC and fully 
consist ent with that of the SMC, i.e., without a 2175 A 
bump (|Fvnbo et al.ll2001a[ ). 

4.2. Gas to dust ratio 

iGordon et al.l (|2003f) suggested that the average ex- 
tinction properties of the SMC, LMC and MW were de- 
scribed by the same underlying extinction law, which 
could be described by the Ry parameter of the MW ex- 
tinction and one additional parameter characterizing the 
strength of the bump and the UV slope. They suggested 
the gas to dust ratio as the second parameter and showed 
that for the LMC and SMC there is a correlation between 
the gas to dust ratio and the UV slope (i.e. 02/ Rv) and 
an anti-correlation between the gas to dust ratio and the 
strength of the bump (i.e. ^bumpMy = ^^/{■y'^Rv)), 
albeit with a large scatter. 

GRBs typically originate in host galaxies with low dust 
to gas ratios (i.e. high gas to dust ratios 'Hiorth et al.J 
2003; Jakobsson et al. 2004c; Vrceswijk et al. 200^ 
Kann et al.l 12006 : iProchaska et al.l 120071 : iSchadv et ail 
2007t ). While the column density, Nh, of the host of 
GRB 070802 is not exceptionally large, it does have a 



TABLE 7 

Extinction and column density of GRBs 



Name 


Av 


logioAf(H I) 


2 


References 


GRB 000301C 


0.09 ± 0.04 


21.2 ±0.5 


2.040 


1 


GRB 000926 


0.18 ±0.06 


21.3 ±0.2 


2.038 


2,3 


GRB 020124 


< 0.2 


21.7 ±0.2 


3.198 


4 


GRB 030323 


< 0.5 


21.90 ±0.07 


3.372 


5 


GRB 030429 


0.34 ±0.04 


21.6 ±0.2 


2.658 


6 


GRB 050401 


0.62 ±0.06 


22.6 ±0.3 


2.899 


7 


GRB 070802 


0.8 - 1.5 


21.5 ±0.2 


2.455 


8 



Note. — The GRBs represented in Figure [TO] as the shaded 
region and GRB 070802. This is not a homogenous sample with 
differing datasets and analysis from the literature. 

Refe rence s. — (1) I Jensen et al.l 1I2OOII') (2) IFvnbo et al.l 
g002al) (3) IFvnbo et al.l COOlal) ~f4') IHiorth et al.l 1I2003|) 

(5) Vrceswiik et al. il20M) (6) IJakobsson et al.l l l2004bli 
(7' . Watson et al.. (.2006i) (8) This paper. 

very large Ay, leading to a low gas to dust. This is 
shown in Fig.[Tni where we have plotted Nn/Ay vs. the 
strength of the bump AbumpMv for the SMC, LMC, 
MW and GRB 070802. Also plotted is the region where 
GRBs for which an extinction curve analysis exists in the 
literature lie. The plot shows that the anti-correlation 
between the bump strength and the gas to dust ratio 
found by'QordoneFin (|200l ) for the SMC and the LMC 
also holds for the GRBs, and GRB 070802 in particular, 
although with large uncertainties in the values. 
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Fig. 10. — Nu/Av vs. Atump/^v- The black filled circles are the averaged values for the SMC (bar and wing (only one line of sight)) 
and LMC (shell and average ) from lGordon et al.l (gOOl). The MW is marked with a star. GRB 070802 is marked by a square, while other 
GRBs from the literature fall into the shaded region bounded by the arrows in the upper left corner. The Nh/Av values for the GRBs 
are calculated from the values quoted in Table [71 Their A^^^^^^/Ay values are estimated to be lower than 0.3 as they are all consistent 
wit h SMC type extinction (i.e. no sign of a bump in the SED). The values for GRB 070802 are calculated directly from the parameters of 
the IFitzpatrick &: Massal II2007I ) parametrization (see Table [Sj. The error bars are dominated by the uncertainty in Ay (and equivalently 
foi' which we have taken Ay = 1.3 ± 0.5 (or Ry = 3.1 ± 1.2). The A-^^^^-p/Ay for the MW is obtained by fitting an FM curve to a 
standard MW extincti on (i.e. Ry = 3.1) and Nu/Ay = 4.93/Rv • lO^^^ cm'^ mag-^ is taken from Diplas & Savage (1994}. We find that 
the trend suggested bv I Gordon et all {2003) of decreasing bump strength with increasing gas-to-dust ratio holds for the MW, GRB 070802 
and the other GRBs from the literature. 

The regions around C I A1560,A1656 are shown in Fig.fTTl 
for GRB 070802 and AO 0235-^164. For comparison we 
also show the same region in the afterglow spectrum of 
the high metallicity burst GRB 000926 which shows no 
significant evidence for the C I absorption features (and 
no bump in its extinction curve). 

It would be of interest to further check whether the 
sightlines to the MW, LMC and SMC for which extinc- 
tion curve analysis exist are consistent with this corre- 
lation. Such a study is beyond the scope of this current 
work, but we note that a quick and incomplete search 
of the literature resulted in two more lines of sight con- 
sistent with this correlation. The first is a detection of 
C I in the line o f sight towards HD 185418 in the MW by 
ISonnentrucker e t al. ( 20 0l). This line of sight is included 
in the sample of .Fitzpatrick fc Massal ()1986i ) and has the 
2175 A bump in its extinction curve. The second is a line 
of sight towards the SMC, which has C I de tected in the 
MW but only as an upper hmit in the SMC (jWeltv et all 
|1997() . Although this sightline does not have an extinc- 
tion curve analysis, it is consistent with our prediction 
that most lines of sight in the SMC should not show sig- 
nificant C I absorption. 

The presence of the C I absorption feature in 
GRB 070802 suggests that the UV radiation field is 
weaker than in typical GRB environments (see § 14. 3p . 



4.3. Detection of C I and the UV radiation field 

GRB 070802 is to our knowledge the first GRB to date 
to show prominent C I in its afterglow spectrum and 
also the first GRB to show the 2175 A feature, which 
is suggestive of a link between these two components. 
Should such a link be established by further observations 
it constitutes a constraint on the environment in which 
the 2175 A feature occurs. The ionization potential of 
neutral carbon is lower than that of hydrogen (11.3 vs. 
13.6 eV) and neutral carbon is therefore only expected 
in regions without intense UV radiation. This would be 
consistent with the general lack of the 2175 A feature in 
GRB host galaxies, as GRBs arise predominantly from 
star-forming regions where the UV radiation is expe cted 
to be strong (jFruchter et al.ll2006l : IChen et al.ll20C)9[ ). 

To test whether this suggested correlation holds in 
general, we have looked at the other three robust de- 
tections of the 2175 A bump in individual systems be- 
yond the local grou p. For the lensing galaxy reported by 
iMotta et al.l (|2002t l and t he intervening absorber towards 
GRB 060418 reported by'Elhson_e£al] (|2006f) . the exist- 
ing data do not cover the wavelength range where one 
would see the C I absorption line if present. However, 
for the intervening damp ed Lyman-g syste r n tow ard 
AO 0235-^164 reported bv I Junkkarinen et al.l (|2004D we 
see tentative evidence for the C I A1656 absorption line. 
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iGordon et al.l (|2003f ) reached the tentative conclusion 
that the shape of the extinction curve is affected by the 
UV flux density in the environment of the dust. In partic- 
ular, a weaker UV flux density is found to correlate with 
the presence of the bump. Continuing our comparison to 
the bump-less GRB 000926 host, we find evidence that 
GRB 070802 does indeed have a weaker UV radiation 
field: i) the GRB 000926 absorption system has stronger 
high ionization lines and much weaker lines from neu- 
tral species (e.g., C I, see Fig.HH) than the GRB 070802 
system, ii) the host galaxy of GRB 000926 appears to 
have a stronger UV flux densit y as illustrated by the very 
strong Lyman-a emission line (|Fvnbo et al.ll2002bD . For 
GRB 070802 we can exclude the presence of such a strong 
Lyman-a emission line (Fig. [T] and Milvang- Jensen et al., 
in preparation) . Thi s supports the conclusion reached by 
IGordon et afl (|2003D . 

It is not immediately clear what the physical signifi- 
cance of the possible correlation of strong C I and the 
2175 A bump is. Given that the bump is generally be- 
lieved to be carried by carbonaceous material, and that 
carbonaceous grain growth a nd formation requires free , 
neutral carbon and molecules (jHenning fc Salamalll998h . 
it would not be surprising to find both observed proper- 
ties in the same environments. The simultaneous pres- 
ence of C I and the 2175 A bump as well as the lower 
overall ionization state of the gas, relatively (though not 
exceptionally) high metallicity, and large dust-to-gas ra- 
tio may be explained in a scenario in which the dust col- 
umn is strongly enriched by the presence of asymptotic 
giant branch (AGB) stars. 

For massive stars to move onto the AGB requires at 
least 60 Myr and typi cally much longer for a large pop- 
ulation (jMaedeil Il992t l. The star- forming environment 
at such an age will be intrinsically relatively benign, 
with a softer UV field, and one in which a large amount 
of dust and molecular and free carbon is pro duced 
(|Andersen et al.ll2003t [Gautschv-Loidl et al.ll2004D . Fur- 
thermore the interstellar medium (ISM) is likely to be 
reasonably metal-rich and dust-rich. These properties 
are in contrast to the normal environments of GRB 
hosts which are typically metal-poor. H owever some 
GRB hosts may be fairly metal-enriched (jFvnbo et al.l 
l2006l ) but still h ave hard radiation fields and young stel- 
lar populations fLe Floc'h et al.|[2003tlChristensen et al.l 
|2004 iProchaska et al.ll2004l) and in particular, low dust- 
to-gas ratios |Fvnbo et al.ll2006t iJakobsson et al.l [20061 
I Watson et al.ll2007h . This is consistent with the above 
scenario since the metal-enrichment timescale could well 
be much shorter than the > 10^ yr required to have 
a reasonable number of AGB stars producing dust 
(jSchneider et al.ll2004[ ). Such a scenario is then also con- 
sistent with the fact that GRB 070802 is the only GRB 
host galaxy so far discovered with a 2175 A bump. It 
should also be noted that the host galaxy of GRB 070802 
is fairly luminous and red for a GRB host (jSavaglio et al.l 
[2009). suggesting that it is a massive, evolv ed system, 
which would be in agreement with the claim of lNoU et al.l 
(|2007t) that the presence of the bump requires an evolved 
population. 

5. CONCLUSIONS 

We have presented VLT observations of the after- 
glow of GRB 070802. In a low resolution spectroscopy 
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Fig. 11. — A comparison of the region around C I in the normal- 
ized spectra of GRB 000926, GRB 070802, and AO 0235-1-164. 
The vertical lines show the positions of the C I lines. The 
two GRB systems have nearly identical Si II line strengths, but 
GRB 000926 has stronger C IV and no significantly detected 
C I. GRB 070802 also has much stronger Al II and Fe II than 
GRB 000926. AO 0235+164, which has a similar H I column den- 
sity as the two GRB sightlines, displays both the 2175 A bump 
and significant C I. 

of the optical afterglow we detect a large number of 
strong metal lines from absorption systems at z = 2.078, 
z = 2.292 and z = 2.4549. The highest redshift system is 
remarkable in showing very strong metal lines, e.g. with 
higher Wr for the Si II lines at 1526 and 1808 A than for 
any other known absorption system we are aware of. We 
also detect strong absorption from C I implying that the 
gas is shielded from strong UV radiation. The spectrum 
shows a red wing of a Lya line from which we derive a H I 
column density of log(Af(H I)) — 21.50±0.20. Imaging of 
the field revealed a fairly bright and red host, detected 
both in R and K bands, suggesting that it is an evolved, 
massive galaxy. 

The spectrum is also remarkable in that the extinction 
curve of the line of sight towards GRB 070802 shows a 
clear signature of the so called Milky Way or 2175 A 
bump. At a redshift of z = 2.45, it is by far the highest 
redshift detection of the 2175 A bump to date. It shows 
that the conditions for the creation and non-destruction 
of the carrier of the bump must already have been in 
place early in the universe. This is the first clear detec- 
tion of the bump in the host of a GRB, with the SMC 
being the typical type of extinction for GRB sightlines. 
This makes GRB 070802 an ideal candidate to study 
the environment needed for the creation and/or non- 
annihilation of bump, by comparing it to other bump-less 
GRBs. 

To accurately derive the properties of the extinc- 
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tion curve in the UV and the bump w e fit it to the 
parametrization of iFitzpatrick fc Massal (|2007l ). The 
bump is found to be centered at Xc ~ 4.6 ±0.1 /im~^ 
(or Ac ~ 2174 ± 50 A), the width of the bump is found 
to be 7 « 1.08 ± 0.05 /im^^ while the 'strength' is 
C3 ~ 2.7 ± 0.1 (taking into account the uncertainty in 
the intrinsic spectral slop e /3). The va l ue of C3 is the 
same as the average that iGordon et all (|2003D find for 
the LMC average sample (C3 = 2.7 ± 0.1) but higher 
than for the LMC2 sample (ca = 1.5 ± 0.1) although it is 
within the scatter. The width 7 is a bit wider compared 
to their 7 = 0.93 ± 0.02 /Ltm"! forethe LMC average and 
7 = 0.95 ± 0.03 for LMC2 (although it is still within the 
scatter of both samples). The amount of extinction is 
Ay « 1.3, but when taking into account 1 a deviations 
in /?, Ay = 0.8-1.5 for the FM and LMC fits. 

In the Local Group, the MW bump is a characteris- 
tic feature of the MW extinction curve. It is also ob- 
served in the LMC, although it is usually weaker and 
followed by a steeper rise in the UV, while it is not ob- 
served in four out of five curves measured for the SMC 
fsee ll.ll for more details and exceptions). Of these three 
'local type extinction laws' we find that the extinction 
of GRB 070802 most closely resembles that of the LMC. 
We find that this result is robust, even taking into ac- 
count a possible contribution to the extinction from the 
strong foreground Mg II absorber. It has been suggested, 
based on the difference in the SMC, LMC and MW, 
that the s trength of the bump correlates with metallicity. 
However. iGordon et al.l (|199 7^ found that although star- 
burst galaxies can have varying mctallicitics, their ex- 
tinction curves all lack the 2175 A bump. By comparing 
GRB 070802 to another high metallicity GRB sighthnc 
which does not show any sign of a bump in its extinction 
curve, we similarly conclude that metallicity is not the 
only driver of the 2175 A bump. 

Another special feature in the spectrum of GRB 070802 
is the detection of a strong C I absorption. This is to our 
knowledge the first GRB spectrum to contain C I ab- 
sorption and we propose that there may be a correlation 
between its detection and the presence of the bump. We 
have checked this suggested correlation for a few other 
lines of sight, and found all of them to be in agreement. 
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Fig. 12.- — A comparison of the SMC and LMC 'typical' extinction curves from'Peil (IT993) and [Gordon et al.l 1 12003) . The solid curves 
are the |Pei] (1992) curves for A^/ = 1, the diamon ds ar e the Gordon ct al. (2003) curves for Ay = 1 for their "LMC2" and "SMG Bar" 
average with error bars, and the dashed curves are iPeil 1(199 2') fits to the Gordon et al.i t2003i') curves (where Ay is varied). The fits are 
done for points in the shaded region which is the same region as for the fits in the paper. The derived Ay differ by 5% for the SMC curve 
and by 2% for the LMC curve which is m uch smaller than the uncertainty in Ay in our fits from the uncertainty in the intrinsic slope /3. In 
addition, the UV slope of the |Pei| l) l992f ) LMC type of extinction is slightly steeper than that of Gordon ct al. (2003), but falls within the 
1 a limit. We therefore conclude that the choice of parametrization of the SMC and LMC type of extinction does not affect our comparison 
of the extinction curve of GRB 070802 to the local type of extinction laws. 

PARAMETRIC EXTINCTION LAWS 
In this Appendix we describe the parametrizations we have used to model the extinction curves. 

The Pei parametrization for the SMC and the LMC 
This parametrization was introduced bv iPeil ()1992[ ) and is given by: 



t (A/A. 



(A./A)". 



(Al) 



where the six terms are all positive and represent different parts of the extinction curve. The parameters can be found 
in ,Peii (,1992. ). The five terms with Ui = 2 are equivalent to Prude profiles with a peak at A^. The only free variable in 
the fit is the overall amount of extinction. The orig inaim (fT99l) pap e r scal es it to but we choose to scale it to 
Ay to be consistent with the other parametrizations. Note that the IPeil (|1992f ) law can also be used to describe Milky 
Way type of ext i nction . 

I Gordon et al.l (|2003f ) present new and updated a verage extinction curves for the SMC and LMC extinction curves. 
Their analysis is based on using the iFitzpatrick fc Massa U9 .9Q.) parametrization which differs from the updated FM 
parametrization we use (see below) in keeping C5 fixed. The iGordon et al] ()2003l ) analysis presents a more nuanced 
picture of the extinction in the SMC and LMC with lines of sight showing different type of extinction (see II. ip . 
However, we find that their "LMC2" avera ge extinctio n curve and their "SMC Bar" extinction curves are very si milar 
to the LMC and SMC extinction curves of IPeil (|1992l . see Fig. IT^ . thus justifying our use of the commonly used IPeil 
([1993) parametrization. 

The CCM parametrization for Milky Way type of extinction 

This parametrization of the Milky Way extinction law was proposed bv iCardelli et al.l (119891 ). It depends on only 
two parameters, E{B — V) — A{B) —A{V) and Ry — A{V)/E{B — V) which is the ratio of total to selective extinction. 
It is given by 



Ax=E{B - V) [Rya{x) + b{x)] 



(A2) 



.A{V) 



ity 



where ^(A) is the total extinction at waveleng th A, a{x) and b{x) are polynomials and x = X ^. The advantage of 
this parametrization over the one of iPeil (|1992f ) for the Milky Way is that it allows for a varying Ry. 



The FM parametrization 
The parametrisation for the UV (i.e. valid for x > 3.7 /im~^) is given by 

Ax=E{B-V) {k{X-V) + Rv) 
= A{V) f±.k{X-V) + l 



(A3) 
(A4) 
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where 



and 



^ Ci+ C2X + C^U{x,XcT^) + Ci{X ~ C<^) X > , 



X2 



D{x, xr_, 7) = — 5-^- — —. (A6) 

The parameters ci and C2 define the Unear component underlying the entire UV range, C3, Xc and 7 give the 2175 
A bump (although its central wavelength is not fixed in the parametrization) and C4 and C5 give a far-UV curvature 
component. The extinction properties in t he infrared and optical are not parametrized in the FM2007 description but 
are derived using spline interpolation (see iFitzpatrick fc Massall2007l for details). As our dataset does not reach into 
these regions in the restframe, we do not constrain the extinction cur ve in this region. Ther e fore, for display purposes, 
we have chosen to set these parameters to 'typical' values found bv IFitzpatrick &: Massal (|2007D . to create 'normal' 
smooth continuation of the curve towards a; = 0. We note that the choice of these parameters does in no way affect 
our fits for the UV p arameters and, vic e versa , that the FM fits do not constrain this part of the curve. 

As explained in IFitzpatrick fc Massal ()2007D . additional useful quantities can be defined using the UV parameters. 
The first one (1) A1250 = 04(8.0 — 05)^ gives the value of the far-UV curvature term at 1250 A and measures the 
strength of the far-UV curvature; (2) Abump = ''■C3/(27) is the area of the bump and (3) E'bump — c^/l'^ is the 
maximum height of the bump above the linear extinction. 



